Length and time scales of a liquid drop impact and penetration into a granular layer 1 

Liquid drop impact and penetration into a granular layer are investigated with diverse 
liquids and granular materials. We use various size of SiC abrasives and glass beads 
as a target granular material. We also employ ethanol and glycerol aqueous solutions 
as well as distilled water to make a liquid drop. The liquid drop impacts the granular 
layer with a low speed (~m/s). The drop deformation and penetration are captured by a 
high speed camera. From the video data, characteristic time scales are measured. Using 
a laser profllometry system, resultant crater morphology and its characteristic length 
scales are measured. Static strength of the granular layer is also measured by the slow 
pillar penetration experiment to quantify the cohesive force effect. We find that the time 
scales are almost independent of impact speed, but they depend on liquid drop viscosity. 
Particularly, the penetration time is proportional to the square root of the liquid drop 
viscosity. Contrastively, the crater radius is independent of the liquid drop viscosity. The 
crater radius is scaled by the same form as the previous paper, (Katsuragi, Phys. Rev. 
Lett. vol. 104, 2010, p. 218001). 
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1. Introduction 



Gra n ular sys t em b e haves sometimes like fluid, b ut sometimes like solid ([Jaeger et all 
(1996h: iDuranl (2000); lAranson fc Tsimringl (2006)). f a collisional impulse is applied 
to a granular matter, the granular matter behaves like fluid at the impact moment. 
After that, it relaxes right away and behaves like solid at a later stage. Thus the gran- 
ular response observation against the impact can prove both of fluid-like and solid-like 
behaviors. In this sense, the granular impact experiment is very useful to study the 
granular matter dynamics. A lot of granular impact experiments using a solid projectile 
have been carried out recentl y . Some of them were interested in the crater morphol- 



ogy flA mato & Williams (1998 ):IUehara et all (j2003MWalsh et al.l ([2Q03): Ambro so et 



2005 b \):\de Vet fc de Bruvnl( 2007 )), and others discussed the dynamics dPica Ciamarra t 
(|2QQ4h : I Ambroso et all (|2005ah:lHou et all (l2005h:lKatsuragi fc Durianl (2007h:lRover"eH 



et 



al 



2QQ2|i ; lNelsQiL-gl-fl£ J d2QQ8[) ; iGoldman fc Umbanhowarl ([ 2 Q08|)). Granu l ar jet faigroddsen fc Shenl 
200lh : lLohse et all (l200ll.l2004h :lBe rgmann et all (l2QQ5h7 Rover et all (l2QQ5h:ICaballero > et all 
([20071 ): Ivon Kann et all (|2010h and eiecta ([Boudet et all (|2QQ6h : iDeboeuf et all (|2009h ) 
creat ed by the impact have been also studied. Numer ical studies have been performed as 
well (|Tsimring fc Volfsonl (|2QQ5h : IWada et all (|2006h ). Thanks to these numerous stud- 
ies, granular impact morphology and dynamics are partly understood. However, all these 
results have concerned the solid projectile impact. For the solid projectile impacts, the 
projectile is much harder than the target granular layer. There was no need to consider 
the projectile deformation effect. Then the next natural question, which is the main focus 
of this paper, is "what happens when the projectile is deformable?" In order to answer 
this question, we have performed the impact experiment of liquid drop and granular 
layer. 

A liquid drop impact to a hard wall or a li q uid pool phenomenon has b e en also studied 

extens i vely (Marmanis fc Thoroddsenl(1996 ): Range fc Feuillebqisl ( 1998 h Throddsen fc Sakakibaral 



dl998h : iBhola fc Chandral (jl999h : iRioboo et all (|2002h : lYarinl (l2QQ6h : IXul feQIOh ). The" 
milk-crown shape by the drop impact is very beautiful and fascinating. Additionally, 
this drop impact spreading and splashing phenomena relate to many industrial tech nolo- 
gies, s uch as ink-jet printing, rapid spray cooling, mixing of solutions, and so on ( Yarinl 
(200^). Although the s l ow pe netration of liquid drops into a powder bed has been re- 
ported ([Hapgood et all ((2002)), the i mpact between liquid drop and granular layer has 



Nefzaoui fc Skurtvsl 



not be en examined until very recently (jDelon et 
(|2010l )). Dalziel and Seaton have studied similar phenomenon of the impact between a 
droplet and a gran ular layer. But their re sult was very limited and focused only on 
resuspension effect (Dalziel & Seaton (l2003h V 



In the previous study ([Katsuragil (j2Q10h ). SiC abrasives and water drops were used 
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to characterize typical crater shapes and to reveal the crater radius scaling. The impact 
experiments with glass beads and some kinds of dr ops have been carried out recently by 
another research group ( Nefzaoui fc Skurtvsl ( 2010h ). They have derived scalings for drop 
deformation and crater diameter. However, these experiments were still limited to reveal 
the overall characteristics of the drop-granular-impact. In order to clarify the details of 
the drop-granular-impact, we have performed experiments with various liquid drops and 
granular layers. 

In laboratory experiments, parameter range we can access is actually very limited. For 
instance, it is hard to make a large liqui d drop since the capilla ry length of usual liquids 
on the earth surface is about mm scale (|de Gennes et all (|2QQ4l )). In our experiment, the 
impact speed made by the free fall is also limited within the order of m/s. While the 
experimental conditions are limited as written above, the drop-granular-impact experi- 
ment can be a key to understand the geological scale events. Actually, low speed impact 
of small solid impactor to very loose gran ular layer resemb l es larg e planetary impacts, in 
terms of the dimensionless similarity law ([Bergmann et all ((2005) ) . In addition, real geo- 
logical i mpactor s in th e space might be fully broken when it hits the other astronomical 
objects (jMeloshl ((1989)). Furthermore, the gravitational force might be smaller than the 
earth surface. Large liquid drops are possible in such environment. Even if the drop size 
is limited in mm scale, the drop-granular-impact might relate to the fossil rain drop s. The 
fossil r ain drops ar e the small pit-like depressions found in fine grain sediment (|Desor 
(|l850h : lMetzl (|l98lh ). The origin of this fossil shape is still a matter of debate. The mor- 
phological study of the drop-granular-impact might give a hint to understand the origin 
of the fossil shape. 

In the next section, the experimental apparatus, materials and procedures are intro- 
duced. In section [3H1 drop deformation and impact dynamics are discussed based on high 
speed video data. In section [321 a static strength measurement method is introduced to 
characterize the strength of the granular layer's surface. We discuss the drop deposition 
behavior observed in small glass beads impacts, using this measurement result. Next, 
the crater morphology is studied based on three dimensional (3D) surface measurement 
data (section l3.3p . After that, typical time scales are measured and discussed using the 
video data (section [4. ip . In particular, we focus on the viscosity dependence of the time 
scales. Finally, characteristic length scales of resultant crater shapes are measured and 
discussed (section 14. 2p . We find the simple approximate relation between crater depth 
and intersectional area. Previously proposed scaling for the crater radius is verified for 
various experimental conditions. Discussion and conclusion are stated in sections [5] and 



2. Experimental 

Experimental apparatus we built is shown in figure HJ A nozzle is mounted on a height 
gauge. The granular target is set below it. Using a syringe pump, liquid is delivered 
slowly (0.2 ml/min.) to the nozzle. A drop grows gradually on the nozzle's tip, and is 
finally pinched off by the gravitational force. Then the drop commences the free fall. We 
vary the nozzle tip inner diameter (0.25, 0.7, 2.4, or 2.5 mm) to control the drop size. 
Radius of the drop Ri made by this procedure ranges from Ri = 1.3 to 2.4 mm. The 
largest value (2.4 mm) is almost the upper limit for a drop under the gravity. We define 
the free fall height h as height difference between the nozzle tip and the surface of the 
granular layer. It is adjusted by the height gauge. The drop elongates vertically up to 
~ 8 mm just before its break. Thus the minimum value of h is chosen as 10 mm, and the 
maximum value of h is 480 mm. Impact speed v is calculated from the free fall relation 
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Figure 1. Schematic image of the experimental apparatus, (a) A liquid drop is released from a 
nozzle tip mounted on a height gauge. The impact is captured by a high speed camera, (b) The 
resultant crater shape is measured by the laser profilometry system composed by a line laser 
sensor and automatic stages. 



v = ^/2gh 1 where g is the gravitational acceleration (q = 9.8 m/s 2 ) . Air drag deceleration 
is negligible in such low speed regime ([Range fc Feuilleboisl (|l998h ). 

Impacts are captured by a high speed camera, TAKEX FC350CL, which takes 210 fps 
video data with 640 x 480 pixels in 8 bit grayscale images. After drop penetration, 
the crater shape is measured by a line laser displacement sensor based profilometry 
system. This system consists of a line laser displacement sensor, KEYENCE LJG030, 
and automatic stages, COMS PM80B-100X (translational motion) and COMS PS60BB- 
360R (rotational motion). An intersectional profile is obtained by the line laser sensor 
at a certain position. Then the sample is shifted using the stage, and the profile is 
obtained at that position. This procedure is repeated. Finally, the intersectional profiles at 
different positions are combined to make a 3D surface depth map. The system has 33 and 
1 /am horizontal and vertical resolutions, respectively. Almost the same system has be en 
applied to measuring the polymer gel surface pattern formation ([Mizoue et al. (l2Q10hV 
Our p revious study of the drop-granular-impact has also used the same system ([Katsu ragi 
(2010)). Data handling and all instruments control are conducted by Lab VIEW software. 

We use SiC abrasives and glass beads for granular targets. The material true densities 
are 3.2 and 2.5 g/cm 3 for SiC and glass beads, respectively. Microscopic photos of SiC 
grains, which reveal that the SiC grains possess nonspherical shape and polydispersity, 
are shown in figures [2^a,b). The central value of the distributed grain diameter D g is (a) 
4 and (b) 50 jim. Glass beads (D g = 50 jim) image is shown in figure Efc). Glass beads 
possess good spherical shape and small polydispersity. We employ D g = 4, 8, 14, 20, and 
50 jim SiC grains (supplied by Showa denko k.k.; JIS R6001 #3000, #1500, #800, #600, 
and #220, respectively) and D g = 5, 18, 50, and 100 jam glass beads (supplied by Potters- 
Ballotini Co., Ltd.). A small vessel (30 mm in diameter and 10 mm thick) is filled with 
grains as a target. Filling and surface flattening of piled grains are performed by hand. 
Grains are piled up on the vessel using a spoon, and then a flat metal plate sweeps the 
pile to make a smooth flat surface. No additional tapping or other compacting process is 
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Figure 2. Photos of grains: SiC abrasives of D g = (a) 4 and (b) 50 /xm, and glass beads of 
D g — (c) 50 /im. SiC grains are nonspherical and possess polydispersity, while glass beads are 
spherical and little polydisperse. 



Material 


D g (/im) 


Pg (g/cm 3 ) 


packing fraction 


SiC abrasives 


4 


0.99 


0.31 


SiC abrasives 


8 


0.99 


0.31 


SiC abrasives 


14 


1.5 


0.44 


SiC abrasives 


20 


1.5 


0.44 


SiC abrasives 


50 


1.5 


0.50 


Glass beads 


5 


1.0 


0.39 


Glass beads 


18 


1.4 


0.55 


Glass beads 


50 


1.5 


0.59 


Glass beads 


100 


1.6 


0.63 



Table 1. Granular layer properties: grain size D g , bulk density p g , and packing fraction of 

target granular layer. 



applied. This simple preparation method makes good reproducibility. The uncertainty of 
the bulk density is a few percent. The impact crater shape also shows good reproducibility. 
Granular properties we use are summarized in table [TJ All packing fraction values have 
less than 0.02 uncertainty. SiC grain's ir regular shape and polydispersity result in a 
very small packing fraction (jRollerl (|l930[ ): ISuzuki et all (|200ll )). Even for spherical glass 
beads, small packing fraction is observed. In such case, grains are too small, and cohesive 
force affects the packing state. This cohesive small glass beads layer is very strong against 
the impact, as described later. In addition, the surface roughness of the granular layer 



brings a lotus effect, which affects the fragility of the impacted drop (jRevssat et al. 



We also varied liquid drop surface tension and viscosity as well as drop size. Ethanol 
and glycerol are utilized to control surface tension and viscosity. Distilled water, 10wt% 
and 100wt% ethanol aqueous solutions, and 60wt%, 88wt% and 100wt% glycerol aque- 
ous solutions are employed to examine the surface tension and viscosity effect. Liquid 
prope rties are listed in table El These values are picke d up from tables in the data 



re pi< 

book rtThe Japan Society of Mechanical Engineer sj (fl983r >L 



3. Results 

3.1. Dynamic morphology 

First, we focus on high speed video data. In figure [U snap shots of various impacts are 
presented. Each row corresponds to different experimental condition. Column indicates 
the impact moment. It is the definition of the origin of time, t = 0. The impact moment 
is identified by the drop deformation in the video data. This identification method is 
not very precise. The frame rate of the video data is not also very fast. The accuracy 
is almost determined by the frame rate as 1/210 s. Column 1-5 represent the successive 
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Liquid Ri (mm) pi (g/cm 3 ) 7 (mN/m) v (10 6 m 2 /s) 



distilled water 1.3 - 2.4 1.0 72 0.89 

ethanol 10wt% 1.9 0.98 47 1.5 

ethanol 100wt% 1.5 0.79 22 1.5 

glycerol 60wt% 2.0 1.2 68 9.6 

glycerol 88wt% 2.0 1.2 65 120 

glycerol 100wt% 1.9 1.3 63 1400 



Table 2. Liquid drop properties (|The Japan Society of Mechanical Engineers] (| 19831 )). Ri, pi, 
7, and v are drop radius, density, surface tension, and kinematic viscosity, respectively. 



photos with 210 fps, right after the impact. Column 6 and 7 show the final penetration 
stage. 

In figures 09(a,b), small water drop impacts are shown. The water drop radius is 
Ri = 1.3 mm, and the free fall height h is (a) 10 mm and (b) 160 mm. One can con- 
firm (a) the sink-type crater and (b) the ring-type crater . They have been reported in 
our previous paper with large drop size (Ri = 2.4 mm) (jKatsurad (|2010h ). The drop 
deformation dynamics are qualitatively similar to the previous result, as well as the final 
crater shapes. This implies that the drop size hardly affects the drop deformation and 
crater morphology, as long as the drop size is in mm order. Much smaller drops of the 
/im or nm size may reveal completely different impact dynamics and crater morphology, 
though. The drop penetration in these processes is really gentle. They never show any 
sudden change. This process is v ery contrastive t o the fakir-impale transition of the drop 
on a rough hydrophobic surface (jRevssat et all ( 20081 )). 

Typical pictures of the glass beads impacts are shown in figures [3fc-e). In all situations, 
inner ring structure can be observed. Even in slow impact case (figure 09(c)), sink- type 
crater cannot be obtained. Another characteristic feature of the glass beads impact is an 
enhancement of the fingering instability. We can confirm the fingering instability and its 
remnant crater rim shape (figures [3j(d,e)). This fingering instability can be often observed 
for D g > 250 fim cases. This result is natural sin ce the surface roughn ess of the granular 
layer should influence the spreading of the drop ([Revssat et all (|2QQ8h ) . 

When the grain size of glass beads is too small, the drop behavior becomes quite 
different. We show the D g = 5 /im glass beads impact in figure 09(f). We are not able 
to see any cratering in this situation. Instead, the water drop spreads on the surface 
of glass beads layer. And it penetrates into the layer very rapidly. Actually, we observe 
this deposition behavior in all free fall heights (h = 10 - 480 mm) for D g — 5 jam glass 
beads. We think this deposition occurs due to the cohesive force among glass beads. 
The principal origin of this cohesive force is capillary bridge effect which is emphasized 
in smaller particles under the laboratory humidity (50 ± 10%) condition. Since glass is 
very hydrophilic material, the capillary cohesion makes a very strong glass beads layer. 
Hydrophilic glass surface also helps rapid penetration. As a consequence, drop deposition 
takes place in this situation. It should be noticed that small SiC grains (D g = 4 /im) 
are able to make various crater shapes, i.e., deposition never happens. This means that 
the SiC grains are hydrophobic compared to glass beads. This consideration is supported 
with the time scale analysis discussed later. Hydrophilicity of target grains affects the 
impact and penetration dynamics. To characterize this effect, we performed another slow 
penetration experiment. A steel pillar was penetrated into a granular layer very slowly. 
The result of this experiment is discussed in the next section. We have actually observed 
the cohesive effect in D g = 18 /im glass beads as well. In that case, very shallow rim and 
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very small central ring or bump were made. These shapes are too thin to be measured 
using our measurement system. Thus we remove these cases (glass beads D g = 5 and 
18 /im) from the characteristic time and length scales analyses below. 

In figure^g), large glass beads (D g = 100 /im) case is displayed. The drop deformation 
and outer rim structure are similar to those of figure [^d). The resultant crater inner 
shape is close to the bump- type (convex one) rather than the ring- type (concave one). 
This t rend (bump -type crater appears in large D g grains impact) is same as the SiC 
result (jKatsurad (|2010h ). 



Next, we show the low surface tension liquid drop cases in figures [3fh,i). A pure ethanol 
(100wt%) drop impacts to (h) SiC of D g = 4 /im and (i) glass beads of D g = 50 /im. As 
can be seen, the drop is no longer stable owing to the fingering instability. Impact inertia 
easily fragments the drop to smaller drops. In such case, resultant crater size evidently 
becomes smaller than stable drop case. Moreover, unstable drop structure is sometimes 
frozen to the irregular shaped bump, as shown in fgiure[3ji). Because we would like to 
focus only on the drop without fragmentation, we remove these data (ethanol 100wt% 
impacts) from the characteristic time and length scales analyses. Ethanol 10wt% impacts 
are more or less stable. We use only the ethanol 10wt% impact data for the analyses 
below. 

Finally, viscous drop impacts are shown in figures [3jj,k). A pure glycerol (100wt%) drop 
impacts to (j) SiC of D g = 4 /im and (k) glass beads of D g = 50 /im. Since the viscosity of 
glycerol is three orders of magnitude greater than water (table [2]), the impact dynamics is 
quite different. The most significant difference is deformability of the drop. The glycerol 
drop deforms itself much less than the water drop (compare, e.g., fgiures [3jd,g, and 
k)). Another characteristic feature of viscous impact is slow penetration. Since the drop 
penetration takes very long time, we cannot obtain the final crater shape using high speed 
video data as shown in figure Efj). We measured the penetration time using another slow 
video camera in such case. And we will discuss the viscosity dependence of the penetration 
time scale in section I47T1 

3.2. Static strength measurement 

As mentioned above, small glass beads layer is very strong due to the capillary bridge 
effect. To characterize this effect, we measured the static strength of the granular layer 
using a universal testing machine (Shimadzu AG100NX). A steel pillar with 5 mm di- 
ameter was penetrated into the granular layer very slowly (0.1 mm/min.). The granular 
samples we used here were completely same as the impact experiment ones. Stroke z and 
stress a applied to the pillar were measured using the testing machine. The obtained 
stroke-stress relations are shown in fiures |4ja,b). In all a data, crossover behavior is 
confirmed. At large z, asymptotic linear relatio n is observed. We t hink this comes from 



the hydrostatic pressure of the granular layer ([Stone et all ([2004 fel lah). This hydrostatic 



regime lasts to z = 3 mm. This means that the bottom wall effect is negligible up to 
this depth. All the crater depth is shallower than this value. On the other hand, steeper 
linear relation appears in very shallow regime (z < 0.5 mm). We assume this is a kind 
of elastic response of the granular layer. The origin of this pseudo elastic property is the 
grains cohesive force due to the capillary bridge effect among grains. This elastic regime 
and its limit stress are actually very small. 

Here, we define the crossover stress value as the breaking stress erg. We first compute 
the linear fitting for large z regime. And, gb is determined by the a value at which 
the data deviate from the fitting. In figure HJc), gb is plotted as a function of granular 
bulk density p g . Each data point is obtained by the average of three independent runs. 
As expected, the gb of D g = 5 /im glass beads layer is much greater than others. 
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Figure 3. Photos of raw data taken by a high speed camera. Each row corresponds to each 
experimental condition as cited on each left bottom part. Specific experimental conditions are 
as follows: (a) SiC of D g = 4 /xm vs water drop of Ri — 1.3 mm with h = 10 mm, (b) same 
as (a) except h = 160 mm, (c) glass beads of D g = 50 /im vs water drop of Ri — 2.2 mm 
with h — 15 mm, (d) same as (c) except h — 120 mm, (e) same as (c) except h = 240 mm, (f) 
glass beads of D g = 5 fim vs water drop of Ri — 2.2 mm with h — 160 mm, (g) glass beads of 
D g = 100 /xm vs water drop of Ri = 2.2 mm with h — 160 mm, (h) SiC of D g = 4 /xm vs pure 
ethanol (100wt%) drop of Ri = 1.5 mm with h — 160 mm, (i) glass beads of D g — 50 /xm vs 
pure ethanol (100wt%) drop of Ri — 1.5 mm with h — 160 mm, (j) SiC of D g = 4 /im vs pure 
glycerol (100wt%) drop of Ri = 1.9 mm with h = 160 mm, and (k) glass beads of D g = 50 /im 
vs pure glycerol (100wt%) drop of Ri = 1.9 mm with h = 160 mm. Column corresponds to the 
impact moment (t = 0). Column 1 to 5 display successive images right after the impact (1/210 
s interval). Column 6 and 7 show later penetration stages. Scale bar (10 mm) is shown at the 
right bottom. As can be seen, various drop deformation and crater morphology are observed. 
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This implies that the small glass beads layer is very strong against the stress. Liquid 
drop impact cannot deform such strong granular layer. That is a reason why the drop 
deposition is observed in small glass beads impacts (figure [3ff)). The gb of D g = 18 pm 
glass beads is about 8 kPa. SiC of D g = 50 pm also shows gb — 8 kPa. These cases 
are marginal ones to make various crater shapes or not, by the drop impact. Therefore, 
we consider gb — 8 kPa is the critical value for our experiment. If the granular surface 
strength gb is greater than 8 kPa, the granular surface is hardly deformed by the liquid 
drop impact. The estimated impact inertia stress by pi = 1 x 10 3 kg/m 3 and v = 3 m/s 
(~ h = 480 mm) is piv 2 = 9 kPa. Where pi is density of the liquid drop. This value 
is close to the measured marginal stress 8 kPa. This critical stress also coincides with 
the single drop stability limit. Around h = 480 mm, water drop fragmentation begins 
to take place, particularly for larger grains cases (D g > 50 pm). We can reduce this 
capillary bridge effect by baking glass beads. But it was hard to completely get rid of 
this effect under our laboratory environment. Except for the cohesive regime (glass beads 
of D g < 18 //in), gb of SiC is greater than that of glass beads (see inset of figure HJ^c)). 
We guess this is due to the grains shape effect. Since SiC grains have irregular shape, the 
grains network connections are stronger than spherical grains. 

Note that this measurement is a kind of quasi static one. We assume this sort of 
measurement is relevant to discuss the granular strength against the impact. In general, 
static and dynamic properties of granular system do not have to agree quantitatively. 
Nonetheless, dynamic characterization using the static property works well in this case. 

3.3. Crater morphology 

After the penetration of liquid drop, final crater shape was measured by the laser pro- 
filometry system. In this section, we discuss the morphology of the crater. 

First, craters by the small water drop vs small SiC grains layer (D g = 4 pm) are shown 
in figure [5l Figures E^a-d) show raw 3D surface depth maps of different free fall heights. 
The depth d is denoted as gray scale. Corresponding radial depth functions d(r) are 
computed using these surface depth maps, and shown in figure Efe). Where r is distance 
from the center of crater. One can confirm the deep sink- type crater at low speed impact. 
And, clear ring-type crater is observed in high spee d impact. This te ndency is completely 
same as large water drop impact case reported by (jKatsuragil ()2QlQh ) . We do not see any 
qualitative difference between large and small water drop impacts, at least under our 
experimental conditions. 

Next, the large SiC grains (D g = 50 pm) and small water drop (Ri = 1.3 mm) impacts 
are displayed in figure [6] The format of the figure is identical to figure [5l Very shallow 
sink- type crater is observed at low speed impact, and bump-t ype crater i s seen at high 
speed impact. This trend is same as large water drop impacts (jKatsuragil (|2010l )) again. 
The only difference between small and large water drop is absolute crater size. For the 
characteristic length scale analysis in section |4]2j we use the crater radius normalized to 
the original drop radius. 

In fiure water drop (Ri = 2.2 mm) and glass beads layer (D g — 50 pm) impacts 
are shown. The crater morphology of glass beads case is a little bit different from SiC 
grains case. Even in low impact speed regime, clear ring-type crater is observed. While 
the crater depth is not very deep (d{r) is always greater than —0.5 mm), inner ring's rim 
height is much higher than SiC impacts. At h = 160 mm, the maximum height is about 
1.5 mm. Moreover, the depth at the center d(0) is greater than 0. Thus this ring- type 
crater has a feature of the bump-type crater as well. The remnant of the drop fingering 
instability can be observed in figure E^c). We observe the asymmetric crater shape at 
the fastest impact case (figure [3(d)). This irregular shape also comes from the fingering 
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Figure 4. Static stress measurement results. A steel pillar (5 mm diameter) was slowly pene- 
trated into a granular layer, and the stroke z and the stress a applied to the pillar were measured, 
(a) Stress vs stroke for various target grains are shown, (b) is the bottom part of (a), (c) Break- 
ing stress as as a function of granular bulk density p g . ctb denotes the crossover between elastic 
and hydrostatic behaviors. The inset shows the blow up of the bottom part. 

instability. Since the glass is very hydrophilic material, water drop can penetrate into 
the glass beads layer faster than recovery of the spherical shape by the surface tension. 
Although the shape is very irregular, corresponding d(r) is computed from the average of 
all radial direction. This procedure enables us to measure the characteristic length scale 
of such irregular crater shape. 

In figure [5J ethanol (10wt% aqueous solution) drop impacts to the glass beads layer 
(D g = 50 fim) are shown. Roughly speaking, the crater structure in figure [8] is similar to 
that of figure E Due to the low surface tension, fingering instability is enhanced and the 
remnant wavy rim structure is also emphasized (figures [5Jc,d)). As mentioned before, 
pure ethanol (100wt%) drop impact causes the drop fragmentation (figures [3jh,i)). Then 
the resultant crater size becomes small. This drop fragmentation effect is behind the scope 
of this paper. We focus only on the drop deformation and the drop-granular-interaction. 
Therefore, we quit measuring the characteristic length scale of the pure ethanol drop 
impact. Namely, we use only the 10wt% ethanol drop impacts data. The surface tension 
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Figure 5. Surface data of SiC grains (D g — 4 /xm) vs small water drop (Ri = 1.3 mm), (a-d) 
Surface height maps are shown as gray scale images, (e) Corresponding radial depth functions 
d(r) are computed from raw height maps (a-d). Qualitative crater structure is similar to the 
large drop case. Sink-type (a) and ring-type (b-d) craters can be observed. 



of 10wt% ethanol solution is 47 mN/m (table [2]). This value is about 65% of the pure 
water's surface tension. This means that it is hard to vary the surface tension widely 
without drop fragmentation. If the surface tension is a little lower than water, drop 
fragmentation relatively easily occurs. 

Finally, we show the glycerol drop impacts to the glass beads (D g = 50 /im), in 
figure [9] In this case, crater shapes are smooth. At low impact speed, almost the fiat 
surface remains. At high impact speed, a simple concave crater appears with a thin 
inner ring. Consequently, corresponding radial depth functions d(r) are calm. This calm 
shape is certainly caused by the viscosity effect. Viscous resistance prevents the drop from 
absorbing surrounding grains. And the drop penetrates into the granular layer very slowly. 
This process makes smooth concave crater shape. To characterize this viscous effect 
quantitatively, we have measured the characteristic time scales. The relation between 
time scales and viscosity is discussed in the next section. 
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Figure 6. Surface data of SiC grains (D g = 50 fim) vs small water drop (Ri = 1.3 mm), (a-d) 
Surface height maps are shown as gray scale images, (e) Corresponding radial depth functions 
d(r) are computed from raw height maps (a-d). Qualitative crater structure is similar to the 
large drop case again. Shallow sink- type (a), ring- type (b,c) and bump- type (d) craters can be 
observed. 



4. Analyses 

4.1. Time scale analysis 

We define two characteristic time scales t& and t p . When the drop hits the granular layer, 
the drop often rebounds and exhibits the dumping oscillation of the drop deformation. % 
is the settling time of this oscillation. After t = the drop penetrates into the granular 
layer. t p characterizes this slow penetration. It is the duration from the impact moment 
to the complete penetration of the drop. In almost all data, tb and t p were measured by 
eye using high speed video data. In some cases, however, t p is too long (t p > 100 s) to 
be measured by high speed video data. In such cases, the impacts were taken also with 
a 30 fps video camera, and t p was measured based on this lower speed video data. 

In figure [10J the measured time scales tb and t p are shown as functions of We. We is 
the Weber number denned as follows, 

We = (4.1) 

7 

Where 7 is the surface tension of liquid drop. In figure [lOj the same experimental con- 
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Figure 7. Surface data of glass beads (D g = 50 /mi) vs water drop (Ri = 2.2 mm), (a-d) Surface 
height maps are shown as gray scale images, (e) Corresponding radial depth functions d(r) are 
computed from raw height maps (a-d). Clear ring- type structure appears in a broad range of 
free fall height (h < 160 mm). At h = 160 mm (c), fingering instability of the drop is recorded 
in the crater rim structure. At the highest free fall height h — 480 mm (d), fingering instability 
affects the shape of the inner bump. 



ditions data (except the free fall height) are connected. Both of t& and t p vary in broad 
ranges; 0.01 < h < 1 and 0.1 < t p < 1000 (s). Nevertheless, it seems that the time 
scales are almost independent of We, in which the principal variable is the impact speed. 
The time scales actually depend on viscosity. Moreover, the time scales depend on the 
granular media. This granular dependence seems to come from wettability difference and 
surface roughness. 

In figure [TTJ we plot average tb and t p vs kinematic viscosity v for SiC (D g = 4 jam) and 
glass beads (D g = 50 jam). Each data point is calculated from the average of different free 
fall height results. As seen in figure [TTT a). t& of SiC is a decreasing function of while 
t\y of glass beads may look almost constant. As shown in figures [3^j,k), viscosity prevents 
the large deformation of liquid drop. Besides, the viscosity represents the dissipation rate 
for the dumping oscillation of drop deformation. As a result, shows decreasing trend 
by viscosity. The gray line in figure [TTla) represents a tentative scaling 4 ~ ^ -1 / 4 . 

In figure fTUb). the penetration time t p is plotted as a function of v. Both of SiC 
and glass beads data show similar increasing trend. The scaling relation we obtain is 
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Figure 8. Surface data of glass beads (D g — 50 /xm) vs 10wt% ethanol solution drop 
(Ri = 1.9 mm), (a-d) Surface height maps are shown as gray scale images, (e) Correspond- 
ing radial depth functions d(r) are computed from raw height maps (a-d). Crater shapes are 
similar to those in figure [3 Drop fingering instability is enhanced due to the low surface tension. 



approximately written as, 

t p ~ v 112 . (4.2) 
According to the Washburn and related recent result, the penetra tion time o f liqui d 



into porous media sh ould be proportional to the viscosity, t p ~ v (jWashburnl (|192lf ): 



Hapgood et all ( 20021 )). However, scaling (j4.2p is quite different from their prediction. 



Although the granular layer is a certain kind of porous media, it is deformable in the 
current experiment. That might be a reason of different scaling exponent. The penetration 
time scale t p strongly depends on granular media as well. The t p of SiC is about one order 
of magnitude larger than that of glass beads. This implies that the SiC is much more 
hydrophobic than glass beads, as mentioned before. 

4.2. Length scale analysis 

Some characteristic length scales are discussed in this section. Although the crater shapes 
are not very simple as shown in figures [3] and [5]|9j we can calculate some relevant length 
scales using d(r). 

First, the vertical length scale is discussed. We define the central depth at r = as d c = 
d(0). The normalized central depth d c jR\ is shown in figurefT^Ta). In order to characterize 
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Figure 9. Surface data of glass beads (D g 



50 /im) 



10 



pure glycerol (100wt%) drop 
e) 



(Ri = 1.9 mm), (a-d) Surface height maps are shown as gray scale images, (e) Correspond- 
ing radial depth functions d(r) are computed from raw height maps (a-d). All crater shapes look 
smoother than other cases. 



this data, we introduce a normalized integration of d{r) as S/Ri 2 = (1/R 2 ) J °° d(r)dr. 
The S/R 2 indicates the compression or dilatation rate of the impacted granular layer. For 
usual solid projectile impact, S/R 2 should be negative. In other words, we expect that 
the impact always compresses the granular layer. In the drop-granular-impact, however, 
S/R 2 can be positive as shown in figure fl2T b). When the drop hits the granular layer, it 
expands and absorbs surrounding grains. Then it penetrates into the granular layer with 
very gentle sedimentation of grains. This process creates a very sparse layer. That is a 
reason why the dilatation is possible by the drop-granular-impact. This dilatation can be 
observed especially for the heavier granular target cases. In such cases, the compression 
by the impact is not so remarkable. 

We find that these two normalized values (d c /Ri and S/R 2 ) behave very similarly. To 
show this similarity, we plot S/R 2 — d c /R\ vs We in figure fl2T c). While the data scatter 
a little, most of all data are around 0. This implies that the crater shape obeys a simple 
relation, 

S~d c R x . (4.3) 

Note that this relation holds both for concave (compressed) and convex (dilated) crater 
shapes. Equation ([4. 3ft means that the rate of compression or dilatation by the impact 
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Figure 10. Time scale data, (a) Drop oscillation settling time tb and (b) drop penetration 
time t p are shown as functions of We. The identical experimental conditions data (except h) 
are connected. In almost all situations, t p is roughly independent of We over two orders of 
magnitude. While tb slightly depends on We, the dependence is very weak. Thus we regard it 
as roughly constant. Each point corresponds to a single impact experiment. 




Figure 11. Averaged time scales (a) h and (b) t p vs kinematic viscosity v. tb for SiC seems to 
decrease as tb ~ Glass beads h does not show obvious trend, particularly in low viscosity 

regime. t p for both cases show the scaling t p ~ v 1 ^ 2 . The t p of glass beads is much shorter than 
that of SiC due to its hydrophilicity. 
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Figure 12. (a) Normalized crater depth at the center d c /Ri and (b) normalized integrated area 
S/Ri 2 vs We. The global trends of d c /Ri and S/Ri 2 are very similar, (c) Subtraction of these 
two quantities is shown vs We. The relation S/Ri 2 — d c /Ri — is roughly satisfied. This relation 
holds both for concave and convex crater shapes. 



can be expressed only by d c and liquid drop radius Ri. Therefore, we regard d c as a 
representative vertical length scale for the drop-granular-impact. Honestly, the d c value 
itself cannot be predicted from the experimental conditions. It does not show simple We 
dependence as shown in figure fT2Ta). The detail analysis to predict d c behavior is an 
open problem. 

Next, the horizontal length scale is discussed. Once we get the d(r) function, it is easy 
to calc ulate the crater r adius R. R is defined by the outer rim radius. In the previous 
report ( Katsura"gil (|2010h ) . we found that the crater radius exhibits the scaling, 



Ri \pi) 



(4.4) 



In soli d projectile granular impacts, similar 1/4 s caling has been reported (jAmato fc Williams 



( 1998[ ): Uehara et al. ( 2QQ3[ ): Walsh et al. (|2QQ3h ). However, the density ratio dependence 
is opposite to the drop-granular-impacts result. On the other h and, We 1 ' 4 " sca l ing fo r 
the water drop defo r mation has been rep orted in the literatures ( Okumura et all ([20031 ): 
IClanet et all (|2004l ) ; iBiance et all (120061)). Thev h ave st udied a water drop impact onto 
a super hydrophobic substrate ( Richard fc Quere ( 2000f )). Since the crater radius seems 
to be determined by the drop deformation degree, the drop deformation scaling is more 
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Figure 13. Scaling plot of the crater radius, R/Ri ~ (p g / ' pi)We 1/A . While this plot includes 
different viscosity results, obtained R data seem to be independent of viscosity. The relevant 
effective parameter to control the crater radius is Ri(p g / p^We 1 ^ 4 . 

plausi ble than t he an alogy of solid projectile impacts, to explain the scaling exponent 
1/4 (|Katsuragil (|2010h ), 

In this study, we varied target grains (size and materials), drop surface tension, and 
drop viscosity. Using these data, we test the universality of the scaling ([4.4)) . In figure HU 
normalized crater radius R/Ri is plotted as a function of (p g / p^We 1 ^. While the data 
are noisy, the global trend follows the scaling (|4.4|) . We can reduce the surface tension only 
to 65% of pure water, using 10wt% ethanol solution. The surface tension variation range 
is actually very limited. The scaling relation ([4.4)) is robust at least in this surface tension 
variation range. Figure [13] includes all different viscosity data. Maybe the glycerol impact 
data are under the scaling line for high We regime, the size dependence on viscosity 
might slightly exist. However, the R/Ri basically looks independent of the viscosity. 
Although the viscosity was varied over three orders of magnitude, the crater radius is 
still independent of viscosity. This result is surprising since the drop deformation depends 
on viscosity (see e.g., figures [3fd,k)). 



5. Discussion 

To vary the crater size, density ratio control is much better than surface tension control. 
Low surface tension drops easily fragment to small drops by the impact. Moreover, the 
surface tension affects the crater radius as 7 -1 / 4 while the density ratio affects as p g /p^ 4 . 
The crater radius is independent of v and is scaled to the impact speed as v 1 ! 2 . Thus v 
or v based control is not so efficient. The density ratio is a quite important parameter 
to control the crater radius. This might mean that the momentum transfer between 
granular layer and projectile is still the most important even for the deformable projectile 
impact. If we assume the momentum transfer as piRfv 2 /di ~ p g RhRi with d\ ~ Rf / 
we obtain the relation Rh/Ri ~ Pg/ ph Where Rh is the horizontal length scale of the 
deformed drop. 

In low viscous drop impact, surface tension governs the drop deformation and naturally 
satisfies ([4.4)) . In high viscous drop impact, we expect that viscosity effect must govern the 
drop deformation. In figures [3jj,k), the drop deformation is suppressed by the viscosity. 
However, the resultant crater radius obeys the scaling ([4.4)) which is independent of the 
viscosity. As long as we observe high speed video data, the deformation rate is quite 
different between water drop and glycerol drop impacts. Nevertheless, the normalized 
crater radius data cannot distinguish them. This means that the data collapse mainly 
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comes from the density ratio effect. The density of viscous drop p\ is relatively large 
(table [2]). Thus the density ratio becomes small and hence drop deformation and R also 
become small. As described above, the density ratio significantly affects the crater radius. 
Therefore, the scaling (|4.4p is recovered even for very viscous drop impacts. 

Contrastively, the characteristic time scales depend on viscosity as £5 ~ z/ -1 / 4 (only 
for SiC) and t p ~ v 1 / 2 (both for SiC and glass beads). They are almost independent of 
the impact speed. The t p scaling is incompatible with traditional Washburn law. A new 
theoretical model might be necessary to solve this problem. 

Since glass is very hydrophilic, the impactor drop can collect a lot of glass beads much 
easier than SiC grains. As a result, the collected many grains produce large inner ring 
or bump for the glass beads impacts. On the other hand, viscosity affects as a resistance 
against the grains absorption. That is why the viscous drop impact cannot make clear 
inner ring or bump. 

Grain size also affects the surface hydrophilicity. And the surface hydrophilicity in- 
fluences the fingering and fragmentation instability of the impacted drop. We have not 
measured the detail of the wettability such as the contact angle. We have used only two 
kinds of grains, SiC or glass. Thus the discussion on hydrophilicity is still qualitative. 
More detail experiments and analyses are required to fully understand the hydrophilicity 
effect. It is the important problem open to future. 

The origin of inner ring f ormat i on m ight be drop interior flow effect, such as the 
capillary flow ( Deegan et al. ( 1997 . 2000f )). May be the fluid flow on the drop surface is 



measureable using absorbed grains as tracers. Bv lClanet et all ([2004t h flow measurement 
has been carried out and the vortical flow has been found for the impact between a liquid 
drop and a hard wall. Flow visualization in the drop is one of the fu ture problems. 

Cha r acterization of drop fingering inst a bility has been also stud ied (Ma rmanis fc Thoroddsen 
(|l996h : iThroddsen fc Sakakibaral (|l998h : iBhola fc Chandral (|l999h ). They have tried to 
identify the boundary between fingering instability and smooth expansion. We observed 
fingering instability in the drop-granular-impact as well. Detail study of this instability 
is also an interesting future problem. 

Some dimensionless numbers made by combinations of Weber and Reynolds numbers 
have been proposed to take into account viscosity effect for the drop impact analysis. Mar- 



manis and Throddsen have derived the impact Reynolds number as / = We^^Re 1 / 2 { Marman is fc Thoroddsen 



1996)). Clanet et al. have proposed Re 1 / 5 scaling for the impact drop deformation (jClanet et al 



20041 )). We tried applying these scaling relations to our data. Unfortunately, they do not 



collapse our data. The simplest one ()4.4j) brings the best result. More detailed investiga- 
tion about the viscosity effect might be required to fully understand the drop-granular- 
impact dynamics. 



6. Conclusion 

We performed the drop-granular-impact experiment using various liquid drops and 
granular layers. We found that t he morpho l ogy o f these impacts is basically similar 
to the previously reported result ([Katsuragil (|2010f )). Grains hydrophilicity affects the 
impact dynamics and crater morphology. Glass beads layer cannot make clear sink-type 
crater. Small glass beads layer has large surface strength due to the capillary bridge 
based cohesive force. This strong granular layer results in the drop deposition rather 
than the cratering. Low surface tension drop impact induces fingering instability and 
fragmentation of the drop. In that case, underlying physical mechanisms are different 
from the single drop impact cratering. Thus we removed the data from quantitative 
analyses. Using the relevant data, we analyzed the characteristic time and length scales. 
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We found that the characteristic time scales are almost independent of impact speed. An 
empirical scaling for the penetration time scale is obtained as, t p ~ v 1 ! 2 . Next, we found 
the simple morphological relation between the impact compression or dilatation and the 
crater depth as S ~ d c Ri. Finally, the previously proposed crater radius scaling R/Ri ~ 
(Pg/ Pi)We x l^ was confirmed in a broad range of the parameter space. Particularly, it is 
surprising that the crater radius is independent of viscosity. 

The author thanks T. Yamaguchi and K. Okumura for useful discussion. This research 
has been supported by the Japanese Ministry of Education, Culture, Sports, Science and 
Technology (MEXT), Grant-in-Aid for Young Scientists, No. 21684021. 

REFERENCES 

Amato, J. C. k Williams, R. E. 1998 Crater formation in the laboratory: An introductory 

experiment in error analysis. Am. J. Phys. 66, 141-143. 
Ambroso, M. A., Kamien, R. D. k Durian, D. J. 2005a Dynamics of shallow impact cra- 

tering. Phys. Rev. E 72 (4), 041305. 
Ambroso, M. A., Santore, C. R., Abate, A. R. k Durian, D. J. 20056 Penetration depth 

for shallow impact cratering. Phys. Rev. E 71 (5), 051305. 
Aranson, I. S. k Tsimring, L. S. 2006 Patterns and collective behavior in granular media: 

Theoretical concepts. Rev. Mod. Phys. 78, 641. 
Bergmann, R., Mikkelsen, R., zeilstra, C, van der Meer, D., Versluis, M., van der 

Weele, K., M, van der Hoef, Kuipers, H. k Lohse, D. 2005 Impact on soft sand: 

Void collapse and jet formation. In Powders and Grains 2005 (ed. R. Garcia-Rojo, H. J. 

Herrmann k S. McNamara), pp. 1211 - 1214. Rotterdam: Balkema. 
Bhola, R. k Chandra, S. 1999 Parameters controlling solidification of molten wax droplets 

falling on a solid surface. J. Materials Sci. 34, 4883-4894. 
Biance, A.-L., Chevy, F., Clanet, C, Lagubeau, G. k Quere, D. 2006 On the elasticity 

of an inertial liquid shock. J. Fluid. Mech. 554, 47-66. 
Boudet, J. F., Amarouchene, Y. k Kellay, H. 2006 Dynamics of impact cratering in 

shallow sand layers. Phys. Rev. Lett. 96 (15), 158001. 
Caballero, Gabriel, Bergmann, Raymond, van der Meer, Devaraj, Prosperetti, An- 
drea k Lohse, Detlef 2007 Role of air in granular jet formation. Phys. Rev. Lett. 99 (1), 

018001. 

Clanet, C, Beguin, C, Richard, D. k Quere, D. 2004 Maximal deformation of an im- 
pacting drop. J. Fluid. Mech. 517, 199-208. 

Dalziel, S. B. k Seaton, M. D. 2003 Resuspension by droplets. In Sedimentation and Sedi- 
ment Transport (ed. A. Gyr k W. Kinzelbach), pp. 101-104. Dordrecht: Kluwer. 

Deboeuf, S., Gondret, P. k Rabaud, M. 2009 Dynamics of grain ejection by sphere impact 
on a granular bed. Phys. Rev. E 79 (4), 041306. 

Deegan, Robert D., Bakajin, Olgica, Dupont, Todd F., Huber, Greg, Nagel, Sid- 
ney R. k Witten, Thomas A. 1997 Capillary flow as the cause of ring stains from dried 
liquid drop. Nature 389, 827-829. 

Deegan, Robert D., Bakajin, Olgica, Dupont, Todd F., Huber, Greg, Nagel, Sid- 
ney R. k Witten, Thomas A. 2000 Contact line deposits in an evaporating drop. Phys. 
Rev. E 62 (1), 756-765. 

Delon, G., Dorbolo, S., Vandewalle, N. k Caps, H. 2009 Drop impact on sand: from 
donut to pie. 62nd Annual Meeting of the APS Division of Fluid Dynamics 54, AH. 00006. 

Desor, E. 1850 On fossil rain drops. Edinburgh New Phil. J. 49, 246 -248. 

Duran, J. 2000 Sands, powders, and grains: An introduction to the physics of granular materi- 
als. New York: Springer. 

de Gennes, P.G., Brochard-Wyart, F. k Quere, D. 2004 Capillarity and Wetting Phe- 
nomena - Drops, Bubbles, Pearls, Waves -. New York: Springer. 

Goldman, D. I. k Umbanhowar, P. 2008 Scaling and dynamics of sphere and disk impact 
into granular media. Phys. Rev. E 77, 021308. 



Length and time scales of a liquid drop impact and penetration into a granular layet21 

Hapgood, Karen P., Litster, James D., Biggs, Simon R. & Howes, Tony 2002 Drop 
penetration into porous powder beds. Journal of Colloid and Interface Science 253 (2), 
353 - 366. 

Hou, M., Peng, Z., Liu, R., Lu, K. & Chan, C. K. 2005 Dynamics of a projectile penetrating 

in granular systems. Phys. Rev. E 72 (6), 062301. 
Jaeger, H. M., Nagel, S. R. & Behringer, R. P. 1996 Granular solids, liquids, and gases. 

Rev. Mod. Phys. 68, 1259. 
von Kann, Stefan, Joubaud, Sylvain, Caballero-Robledo, Gabriel A., Lohse, Detlef 

& van der Meer, Devaraj 2010 Effect of finite container size on granular jet formation. 

Phys. Rev. E 81 (4), 041306. 
Katsuragi, Hiroaki 2010 Morphology scaling of drop impact onto a granular layer. Phys. Rev. 

Lett. 104 (21), 218001. 

Katsuragi, H. & Durian, D.J. 2007 Unified force law for granular impact cratering. Nature 
Phys. 3, 420-423. 

Lohse, D., Bergmann, R., Mikkelsen, R., Zeilstra, C., van der Meer, D., Versluis, M., 

van der Weele, K., van der Hoef, M. & Kuipers, H. 2004 Impact on soft sandivoid 

collapse and jet formation. Phys. Rev. Lett. 93, 198003. 
Lohse, D., Rauhe, R., Bergmann, R. & van der Meer, D. 2001 Creating a dry variety of 

quicksand. Nature 432, 689-690. 
Marmanis, H. & Thoroddsen, S. T. 1996 Scaling of the fingering pattern of an impacting 

drop. Phys. Fluids 8, 1344-1346. 
Melosh, H. J. 1989 Impact Cratering. New York: Oxford University Press. 
Metz, R. 1981 Why not raindrop impressions? J. Sediment. Petrol. 51, 265 -268. 
Mizoue, T., Aoki, Y., Tokita, M., Honjo, H., Barraza, H. J. & Katsuragi, H. 2010 

Control of polymer gel surface pattern formation and its three dimensional measuremnt 

method. J. Polym. Eng. 30, 523 - 534. 
Nefzaoui, E. & Skurtys, O.R 2010 Impact of a liquid drop on a granular- mediu: inertia, 

viscosity and surface tension effects on the drop deformation. ArXiv:1009.1806vl. 
Nelson, E. L., Katsuragi, H., Mayor, P. & Durian, D. J. 2008 Projectile interactions in 

granular impact cratering. Phys. Rev. Lett. 101, 068001. 
Okumura, K., Chevy, F., Richard, D., Quere, D. & Clanet, C. 2003 Water spring: A 

model for bouncing drops. Europhys. Lett. 62, 237-243. 
Pica Ciamarra, Massimo, Lara, Antonio H., Lee, Andrew T., Goldman, Daniel I., 

Vishik, Inn A & Swinney, Harry L. 2004 Dynamics of drag and force distributions for 

projectile impact in a granular medium. Phys. Rev. Lett. 92 (19), 194301. 
Range, K. & Feuillebois, F. 1998 Influence of surface roughness on liquid drop impact. J. 

Colloid Interface Sci. 203, 16-30. 
Reyssat, M., Yeomans, J. M. & Quere, D. 2008 Impalement of fakir drops. EPL 81, 26006. 
Richard, D. & Quere, D. 2000 Bouncing water drops. Europhys. Lett. 50, 769-775. 
Rioboo, R., Marengo, M. & Trope A, C. 2002 Time evolution of liquid drop impact onto 

solid, dry surfaces. Exp. Fluids 33, 112-124. 
Roller, P. S. 1930 The bulking properties of microscopic particles. Ind. Eng. Chem. 22, 1206- 

1208. 

Royer, John R., Corwin, Eric I., Eng, Peter J. k, Jaeger, Heinrich M. 2007 Gas- 
mediated impact dynamics in fine-grained granular materials. Phys. Rev. Lett. 99 (3), 
038003. 

Royer, J. R., Corwin, E. I., Flior, A., Cordero, M.-L., Rivers, M. L., Eng, P. J. k, 

Jaeger, H. M. 2005 Formation of granular jets observed by high-speed x-ray radiography. 

Nature Phys. 1, 164-167. 
Stone, M. B., Barry, R., Bernstein, D. P., Pelc, M. D., Tsui, Y. K. & Schiffer, P. 

2004a Local jamming via penetration of a granular medium. Phys. Rev. E 70, 041301. 
Stone, M. B., Bernstein, D. P., Barry, R., Pelc, M. D., Tsui, Y. K. & Schiffer, P. 

20046 Getting to the bottom of a granular medium. Nature 427, 503 -504. 
Suzuki, M., Sato, H., Hasegawa, M. & Hirota, M. 2001 Effect of size distribution on tapping 

properties of fine powder. Powder Tech. 118, 53-57. 
The Japan Society of Mechanical Engineers, ed. 1983 JSME Data Book: Thermophysical 

Properties of Fluids. Tokyo: JSME, in Japanese. 



22 



Hiroaki Katsuragi 



Thoroddsen, S. T. & Shen, A. Q. 2001 Granular jets. Phys. Fluids 13, 4-6. 

Throddsen, S. T. & Sakakibara, J. 1998 Evolution of the fingering pattern of an impact 

drop. Phys. Fluids 10, 1359-1374. 
Tsimring, L. S. k, Volfson, D. 2005 Modeling of impact cratering in granular media. In 

Powders and Grains 2005 (ed. R. Garcia- Rojo, H. J. Herrmann & S. McNamara), pp. 1215 

- 1218. Rotterdam: Balkema. 
Uehara, J. S., Ambrosso, M. A., Ojha, R. P. & Durian, D. J. 2003 Low-speed impact 

craters in loose granular media. Phys. Rev. Lett. 90, 194301. 
de Vet, S. J. & de Bruyn, J. R. 2007 Shape of impact craters in granular media. Phys. Rev. 

E 76, 041306. 

Wada, Koji, Senshu, Hiroki k, Matsui, Takafumi 2006 Numerical simulation of impact 

cratering on granular material. Icarus 180 (2), 528 - 545. 
Walsh, A. M., Holloway, K. E., Habdas, P. & de Bruyn, J. R. 2003 Morphology and 

scaling of impact craters in granular media. Phys. Rev. Lett. 91, 104301. 
Washburn, Edward W. 1921 The dynamics of capillary flow. Phys. Rev. 17 (3), 273-283. 
Xu, Lei 2010 Instability development of a viscous liquid drop impacting a smooth substrate. 

Phys. Rev. E 82 (2), 025303. 
Yarin, A. L. 2006 Drop impact dynamics: splashing, spreading, receding, bouncing... Ann. Rev. 

Fluid Mech. 38, 159-192. 



